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THEORETICAL STUDIES OF TONE NOISE FROM A FAN ROTOR 

By G.V.R. Rao, W.T. Chu, and R.V. Digumarthi 
Rao and Associates, Inc. 


SUMMARY 

Noise sources in a fan rotor due to inflow distortion 
and turbulence are examined. Analytical studies show that 
high levels of tone noise at blade passing frequency and 
its higher harmonics can result from relatively small dis- 
tortions of the inflow to the rotor. It is shown that these 
discrete tones can also originate from blade loading fluctu- 
ations caused by inflow turbulence, provided the longitudinal 
velocity correlation length scale is sufficiently large. 
Computations carried out on a small scale subsonic rotor 
indicate that acoustic radiation from blade thickness effect 
and fluctuating Reynolds stresses is small compared to that 
from blade loading. Theoretical predictions of on-axis noise 
are in reasonable agreement with acoustic measurements taken 
on a small scale rotor. 
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1-0 INTRODUCTION 

The various mechanisms by- which -noise is generated in a 
lift fan have been the subject of an analytical study during 
the past two years, at Rao and Associates, Inc., under a con- 
tract with National Aeronautics and Space Administration. 

"Lift fans operate essentially on the same principles .as a. 
single stage compressor or a turbo-fan, and the acoustic 
studies on such machines , carried out by several investiga- 
tors over the past decade, guided us in our investigations. 
Noise measurements on a small scale f an ^ 1 ^ and on the various 
OGV configurations of the 36 inch diameter GE LF336 lift fan 
indicated high level discrete tones, contrary to theoretical 
predictions of blade-vane interaction effects. . Consequently, 
the noise sources associated with the rotor itself are inves- 
tigated further and some preliminary results were presented 
• . ■ ( 3 ) 

m our interim report on the contract. 

The noise sources, whether they be from forces or flow - 
interactions can be classified into two general catagories: 
ordered and random. In the first catagory, the sources at any 
one point in the flow region have a deterministic relation , 
with those at any other point. When the sources are random, 
as in the second catagory, we can only mention the spectral 
distribution of the source strength and correlation over lim- 
ited separation distances. Ordered sources, except in a spec- 
ific type of inflow distortion, radiate zero tone-noise to 



points on rotor axis, contrary to the high levels observed in 
test data, and therefore emphasis is given in this report to 
random noise sources and radiation to on-axis points. 

Since the noise generated in a fan is governed by the 
velocities occurring through the fan, our report begins with 
a description of flow field as given in section 2. The con- 
sequences of the flow, such as forces on blades, effect of 
their thickness and the occurrence of fluctuating Reynolds 
stresses are also discussed in this section. The presence 
of a rotating blade, i.e. its loading and thickness, are re- 
placed by suitably phased elemental dipoles and monopoles 
around the rotor circumference, and sound radiation from all 
such sources in the fan annulus is examined in the following 
two sections. The influence of inflow distortion is discussed 
in section 3 and the effect of turbulence in the inflow to the 
rotor is treated in section 4. Acoustic radiation from fluctu- 
ating Reynolds stresses is evaluated in these sections using 
the space-time correlation functions of the velocity fluctua- 
tions. Our primary effort is on estimation of noise sources 
occurring in a fan-rotor, and the effect of duct cut-off has 
been ignored in the preliminary noise calculations carried 
out on the small scale fan and presented here. It is ob- 
served that the random loading on the rotor blades caused by 
inflow turbulence, with sufficiently large length scale, can 
give rise to discrete tone noise at blade passing frequency 
and its higher harmonics. The configuration of the small 
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scale rotor used in our acoustic computations is described in 
reference 1, and the pertinent design parameters are given 
in Fig. 1.1. 

Inflow turbulence and on-axis noise measurements obtained 
by NASA Ames Research Center on the small scale rotor are pre- 
sented in section 5. Acoustic computations, at the flow para- 
meters corresponding to the test conditions, show reasonable 
agreement with measured noise levels. Conclusions and recom- 
mendations for future investigations are given in the final 
section of this report. 

The theoretical investigations and acoustic computations 
were carried out at Rao and Associates, Inc., for NASA Ames 
Research Center under contract NAS2-6401. The helpful sug- 
gestions by Mr. B.K. Hodder and Mr. D.H. Hickey, the technical 

managers on this contract, are greatly appreciated. 

) j’ 
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2.0 


NOISE SOURCES IN A ROTOR 


A comprehensive theory of generation and propagation of 

( 4 ) 

aerodynamic sound is given by Lighthill in his 1961 Bakerian 
Lecture. Following his approach, the rotor noise problem can 
be investigated by considering the role of the following: 

(a) Aerodynamic forces on the blades giving rise to acoustic 
dipoles. 

(b) Effect of blade thickness giving rise to acoustic monopoles. 

(c) Fluctuating Reynolds stresses, and the resulting quadrupoles. 

The monopole and dipole sources are located at the blades, 
whereas the quadrupole sources are present in the entire flow 
region exterior to the blades. All the above noise sources 
are directly related to the velocities occurring in the region 
of the rotor and hence a description of the flow field is 
given in this section, followed by a general discussion of the 

r 

sources and their locations . The coordinate system employed 
in describing the flow field and blade locations is shown in 
Fig. 2.1. 


2.1 Description of Flow Field 

The velocity field in a fan, including the effect of 
moving rotor blades and stationary vanes or struts has been 
described in detail in reference 3. Since the present inves- 
tigation is limited to noise from the rotor itself, we shall 
limit our attention only to the effect of rotor blades. The 
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flow conditions within the wake of each blade have consider- 
able random components, unrelated to the conditions occurring 
in the wakes of the neighboring blades , as pointed out in 
reference 5, and such velocity fluctuations by themselves do 
not contribute to discrete tone noise. Consequently, we shall 
exclude the region downstream of the rotor and describe the 
velocity vector in the flow region as 

V = U + IJ + u (2.1) 

- P 

The inflow to the fan, denoted by IJ, may be both unsteady and 
nonuniform. The velocity resulting from the potential flow 
around the rotating blades is represented by U . The last 
term u denotes the turbulent fluctuations of velocity in the 
inflow. We note that each term, except u, in the above equa- 
tion is well correlated, i.e. the value at any one point is 
definable in terms of its value at any other point. On the 
other hand, the turbulent fluctuations u may have correlation 
only in a small region. 

In the following discussion, we will be describing the 
velocities in terms of their components in the axial, radial 
and tangential directions indicated in Fig. 2.1, by using sub- 
scripts 1, 2, and 3 respectively. Subscripts r and i will 
be used to denote respectively the components along and normal 
to the mean relative velocity met by the rotor blade. 
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2.1.1 Inflow to the fan 


Due to disturbances upstream of a fan engine, the inflow 
to the rotor can exhibit spatial nonuniformity as well as 
fluctuations of such distorted flow. At any given location 
in the fan annulus, the inflow velocity U can be described by 
its components in the axial, radial and tangential directions. 
However, the fluctuations of lift loading on the blades is 
primarily dependent upon the component normal to the relative 
velocity. Denoting this component of the inflow by U x , we 
can express its periodic fluctuations as 

U j (r ,e ,t ) = £ U. A (r,0)*e n n (2.2) 

h 3 

We note that the left hand side is a real function in the above 
representation. Here and in the rest of this report, when no 
limits are shown, the summation is over all integer values 
from - 00 to +°° including zero. The steady state value of U x 
is the term for h = 0 . The variation 1 in the 0 direction 
can also be represented by a similar complex Fourier series 
to allow the following description 


i(gO+<p ) -i (u>. t+<J>, ) 
u i <r,6 ,t) -EJ v h Ve » -e h h 

g h y 




Uge- ah t+* gh 


(2.3) 


Above description of velocity allows different functions of t 
to represent the temporal fluctuations at different positions 
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in the annulus. However, the frequency 10 ^ of the ghth com- 
ponent is independent of the coordinates r,0. This require- 
ment is consistent with ignoring the random disturbance in 
the present treatment. The amplitudes U ^ ^ and phase angles 
<f> ^ can be determined from Fourier analysis of flow field sur- 
vey at the fan inlet. The lowest frequency employed in 
t he a bove equation corresponds to the longest period one finds 
from the fluctuations at various locations in the annulus. We 
note that Eq. (2.3) represents uniform axial inflow conditions 
when 0 )^ and g are both zero. 

In considering the effect of blade thickness, by suitable 
source-sink distribution, we require the knowledge of the vel- 
ocity component along the blade chord, which can be approxi- 
mated by the component taken along the mean relative vel- 

ocity direction. Similar to Eq. (2.3) we can deduce the 
expression 


U (r,0,t) = £ X) 

g h 


u 


gh, r 


e i( g Q ~^h t + % h ) 


(2.4) 


2.1.2 Potential flow around blades 

Let us consider an elemental annulus 2Trrdr of the fan 
containing B number of equally spaced identical blades. The 
potential flow around these blades can be obtained in terms of 
two-dimensional flow past a cascade of infinite number of air- 
foils. Each airfoil in the cascade can be replaced by suitable 
distribution of vortices, sources and sinks to evaluate the 
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flow field due to lift on the blades and the effect of their 
thickness. Such analysis was presented in reference 3, and 
it was found that a simplified description of the flow field 
upstream of the leading edges is possible. The perturbations 
from the mean flow are relatively small in the region down- 
stream of the rotor leading edge and hence will be ignored in 
the noise calculation presented in this report. 

Based on reference 3, we can write the following expres- 
sion for the periodic fluctuations of the axial component of 
the rotor-induced velocity in the region upstream of rotor 
leading edge . 


U = £ U 

Pi „ P 1 j n 


(2.5) 


„ " | 2 nng /d| 

where U „ = U *e *exp{i <j> +inB( 0 -ftt)} 

pi 3 n pi,n T Pi,n 

U j <p = the amplitude and phase angle respectively, of 

Pi* Pi.> n the component at blade leading edge 

r) = axial distance upstream from rotor leading edge 

and d = blade spacing. 

The amplitude n and the phase angle ^ n depend upon the 

blade loading, rotor geometry, and blade thickness distribu- 
tion as discussed in reference 3. 

The above simple representation for n is possible 

upstream of rotor leading edge, since phase angle <j>^ n does 
not depend upon g in this region. Computations at pitch 
radius, on the small scale fan, for the corresponding tangen- 
tial component indicate 
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and 


U = U 

p 3 j n Pi,n 


b = 1 + d> 

2 v pi > n 


Consequently, we can write 


iir/2 


U = U *e 

P 3 i n P 1 J n 


( 2 . 6 ) 


The radial components of will be ignored, since the poten- 


tial flow field past the rotor blades is evaluated under two- 
dimensional flow assumptions. 


2.1.3 Inflow turbulence 

The velocity perturbations discussed in preceeding sub- 
sections are well-ordered in the sense that the velocity at 
one point of the flow region has a definite space-time relation 
with the velocity at any other point of the flow region. On 
the other hand, the random velocity fluctuations present in 
the flow entering the fan can be referred to only in terms of 
statistical quantities such as mean square values, probability 
density functions, correlation functions and power spectral 
density functions. In the study of fan noise problems, another 
important quantity is the correlation length l, which is 
defined as the area under the non-dimensional correlation 
function and thus represents some average size of the eddy 
within which two points in the flow field will have velocity 
fluctuations that can be considered as coherent or related. 

The spectral distribution and the correlation length 
are not independent of each other. This is obvious for the 



reason that if low frequency fluctuations are dominant, the 
extent to which two points in the flow field are related will 
be larger. The way they are related depends on the particular 
model chosen to describe the turbulent field, as discussed in 
reference 3. 

Let us assume that the turbulent fluctuations in the 
inflow are homogeneous, isotropic and are characterized by a 
longitudinal velocity correlation function of the type 
exp (- | £ |/£) , where £ is the separation distance. Choosing 
"frozen-convected turbulence" hypothesis, the spectra of the 
components u and u^ can be related to length scale i and the 
convection velocity. From geometrical relation one can derive 
the spectra for u r and u I , as functions of spectral inten- 
sities of u^ and u^ and the blade stagger angle A. In our 
present investigations, we can use the mean axial inflow 
velocity U as a close approximation to the convection velo- 
city. Furthermore, the blade stagger angle being close to 
tt/2 in general, we shall assume that u x = u and u - u , 
in describing the spectral intensities of u x and u r - Con- 
sequently we used 

$ (u) = u 2 -4— {1 + (iui/U ) 2 }~ 1 (2.7) 

11 , 'flu 3 

^ 3 

and 

<J> (w) = u 2 ~ r { 1+3 ( £u)/U ) 2 } { 1 + (£w/U ) 2 } -2 (2.8) 

11 Zl\U 3 cl 

X 3 

in investigating the blade loading and thickness effects 
caused by inflow turbulence. 
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For the quadrupole calculations , as carried out in sub- 
section 4.3, the space-time correlation function of u^ is re- 
quired. Without assuming isotropy, the space-time correlation 

function of u for "frozen-pattern" convected with velocity U 
1 a 

can be written as 

R (?,t) = u 2 *exp 
u ~ 

1 

where the subscripts on i and £ denote the length scales and 
separation distances in the respective directions. If one is 
interested in that part of the noise generated by the true 
temporal fluctuation of the turbulent field, an additional 
time function has to be incorporated as discussed in refer- 
ence 3 . 

2.2 Rotor Blade Loading 

Each blade element, as it travels around the circum- 
ference with a tangential velocity , experiences periodic 
and random fluctuations in loading caused by the inflow dis- 
tortions and turbulence discussed in the preceeding subsection. 
Typical of the airfoil profiles employed for the blades, the 
drag can be ignored in comparison to the lift force. Let us 
consider a blade element of span dr, at radius r, and experi- 
encing a lift force Ldr, where 

L -? p o V r C L= (2 ' 10) 

The lift L per unit span of the blade depends upon the rela- 
tive velocity V , chord c, and lift coefficient C T . 

3 . 


£, 

T, 


(2.9) 
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Let us assume that all the B number of blade elements 
are equally spaced around the circumference 2-rrr and have the 
same cross-section and stagger angle. Let subscript j identify 
the j'th blade, as j takes on integer values increasing in the 
direction of rotor rotation from 1 to B. The coordinate system 
employed in our analysis and the location of the jth blade are 
typically sketched in Fig. 2.1. For the sake of convenience, 
the reference meridional plane 8=0 is chosen so as to contain 
the point, at which the far field sound pressure is evaluated 
in later sections. Even though all the blade elements are 
identical, they may not simultaneously experience the same lift 
force as discussed in later sections 3 and 4: To relate the 

point of application of load on each blade element toT.ts 
position in the annulus, we assumed that the locus of mid-chord 
locations along the blade span is a radial line. Since the 
center of gravity of a blade section is not too far from its 
mid-chord location, such an assumption is a reasonable approxi- 
mation . 

The force on the jth blade can be evaluated according to 
Eq. (2.10) with the appropriate values of V and C T as experi- 

enced by the jth blade. Considering the force L. as concen- 

3 

trated at a point on the j'th blade its coordinates can be 
derived from Fig. 2.1 as 

x . = -XL ' cos! 

3 

r . = r 


and 


8 = 9, + (j - 1) ^ (2.11) 

v D L 
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where X = blade stagger angle, 

t, "*= center of pressure measured from midchord along 
the chord, positive towards the trailing edge, 

and 0 j = location of midchord of "first" blade at time t. 

Furthermore, if the blade with its midchord point cros- 
sing the meridional plane 0 = 0 at time t = 0 is denoted as 
the "first" blade, we can write 

8 7 = «t + U - 1) (2.12) 

3 u r 

where ft = rotor circular frequency. 

The dependence of blade loading on the inflow conditions 
and the resulting sound radiation is discussed in sections 3 
and 4 of this report. 

2.3 Volume Sources Due to Blade Thickness 

The motion of the rotor blade, due to its thickness, 
causes air to be pushed out and drawn back respectively near 
the leading and trailing edge regions. Velocity perturba- 
tions normal to the rotor disc, caused by the finite thickness 

of blades rotating in uniform steady flow, were considered 

. ( 6 ) . , 
by Deming to evaluate the resulting acoustic radiation. 

In the following, we shall consider the general case of the 
blades immersed in non-uniform unsteady flow and derive ex- 
pressions for the fluctuating volume sources due to the effect 
of blade thickness. Similar to the discussion presented in 
the previous subsection 2.2, let us consider an element of 
span dr at radius r on the jth blade. The two-dimensional 
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flow past the blade profile can be approximated. by the follow- 
ing distribution of source-sink strength along the chord. 


= V r .§ (2.13) 

where £ = distance measured from midchord, along the chord, 
positive towards trailing edge, 

V = relative velocity, considered parallel to blade 
chord , 

and b = blade thickness. 

The subscript z,j on the left hand side of above equation de- 
notes the function as dependent on £ measured on the jth 
blade . 

Taking into account the span dr, the j'th blade element 
can be considered as volume sources 


g . = V • dc dr 
4 CJ r dc ^ 


(2.14) 


located at every point -c/2 £ C c/2 over the chord of this 
jth blade. The coordinates of this point £ along the chord of 
the jth blade element can be written, similar to Eq . (2.11), 

as 

x„ . = -£ cosX 
?<7 

r r . = r 

and 0_ . = fit + (j' - 1)2 tt/B - (c sinA)/r (2.15) 

Cj 

Sound radiation from the volume sources described above is 
investigated in sections 3 and 4 of this report. 
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2.4 "Fluctuating Reynolds Stresses" in the Flow Region 

The manner in which fluctuating flows can generate noise 

(7 ) 

as quadrupoles sources was originally formulated by Lighthill 
For subsonic flow, the dominant source term is the "fluctuating 
Reynolds stresses" 


T. . = p V. V . 
il o x.j 


(2.16) 


where the subscripts i,j indicate the components in the three 
dimensional cartesian coordinate directions 1, 2, and 3. To 
obtain the fluctuating Reynolds stresses T^ ^ in the flow region 
through the rotor, we substitute the instantaneous velocity 
vector V as defined by Eq. (2.1). The resulting acoustic 
radiation is discussed in the appendix, and application to 
rotor noise under various flow conditions is presented in the 
following sections 3 and 4. 
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3 -0 ROTOR NOISE GENERATED BY INFLOW DISTORTIONS 

A general description of inflow distortion is given in 
subsection 2.1.1 in terms of its ghth components. Without the 
knowledge of amplitude, frequency and phase angle of each ghth 
component, it is difficult to evaluate the noise generated by 
inflow distortion. However, we can proceed to examine the 
sound radiation resulting from blade loading, thickness effect 
and fluctuating Reynolds stresses in terms of a ghth component 
of the inflow distortion. 

3.1 Acoustic Radiation from Periodic Blade Loading 

Let us consider an element of span dr located at radius 
r on the jth blade and subjected to ghth. component of inflow 
distortion described in subsection 2.1.1. Since the velocity 
fluctuations normal to the blade chord have the dominant influ- 
ence on the lift fluctuations, it is sufficient to examine the 
influence of U ^ ^ described by Eq. (2.3). The lift per unit 
span of the blade given by Eq. (2.10) will be then periodic 
due to the periodic fluctuations of the incidence angle met by 

the blade, and can be related to its amplitude and the Sears' 

. ( 8 ) 

lift response function * Denoting the blade that crosses the 
reference meridional plane at time t = 0, as the "first" blade, 
we can express the fluctuating lift on the j th blade element as 

A 

L gh,j ^^-expHt^-^t + i^T.+iy (3.1) 
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where L gh = p </ V r c ‘V,i -l 5(Y) 


and 


S(y) = Sears' lift response function 
y = reduced frequency = (ok - gSl) 


2V 


In the above representation the phase angle <J> ^ includes 
the phase angle of the gTith component of U ^ L as in Eq. (2.3) 
and also the phase angle in the Sears' function £( y) . We 
note that L and <p ^ are not subscripted by j since all the 
blade elements at radius r and the flow conditions met by 
them are assumed identical. However, the lift fluctuations 

on the jth blade element occur at time t . earlier than those 

0 

on the first blade, where 


Denoting the point of application of the force L 


(3.2) 


gh,3 


by 


t'fri measured along the blade chord from its midchord loca- 
tion, we can use Eqs . (2.11) and (2.12) to obtain its coordin- 

ates, as 

x . = - c cos X 
3 *gh 


and 


r . = r 

J 

9 . = ftt + (<7 - D2 tt/B - 
0 ^ 


Cfr sinA 


(3.3) 


We note that the location l' , on the blade chord is not. sub- 

gh 

scripted by j , since all the blades are considered identical. 

Considering the lift force as concentrated at the point 
defined by its coordinates given in Eq. (3.3), we can employ 
complex Fourier series representation to obtain the force over 
an elemental area rd8 dr at (-; ' ^ cos X , r, 9 ) as ~ 
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h 

6 L , . = -2— Y' exp{-i (w, - gQ + mil) t + im0 + id> , 

gh,g 2-rr ^ h y 

+ i sinA ) /r - i (m - g) (j - 1) 2 tt/b} dr d0 (3.4) 

y 

We note that 5L , . denotes the contribution due to lift on 

the j th blade element and its motion in a circle with rota- 
tional frequency fl. Sound pressure at the far field point 
due to the fluctuating force described above can be obtained as 


6p 


gh t j 


4irSa 


{ (-sinA cos'l* + cos A sin0 . sin'H-SL , .} 

at «? gh,o 


(3.5) 


where the square bracket indicates that the derivative is 
evaluated at retarded time t - S/a Qj and S is the distance of 
the field point from source located at coordinates (- 5 ^ cosA^r.,0) 
Substitution of Eq. (3.4) into the above is carried out in sec- 
tion A2 of the appendix, resulting in Eq. (A2.1) for y 

Considering the loading on all the blade elements contained in 
the fan annulus, we obtain the far field acoustic pressure as 


2 7T 


J gh 


-f fk 


Sp 


gh, j 


(3.6) 


The superscript d and the subscript gh on the left hand side 
of Eq. (3.5) and (3.6) denote that the expression is for the 
contribution from dipole radiation from the effect of the ghtb. 
component of the inflow distortion. After completing the in- 
tegration over 0 and summation over j in the above equation, 
it is shown in section A2 of the appendix that sound radiation 
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at blade passing frequency occurs only when = 0, i.e. 
for steady state inflow distortion. Integration over r can 
be simplified in terms of the value of the integrand at a 
representative radius r g and it is shown in section A. 2 of 
the appendix that the rms value of the nth harmonic of the 
far field sound pressure given by Eq. (3.6) reduces to 


~d 

p 

gn 


K „ 

nB A 

^ «»b f. v t .-r h ) (3.7) 


where K = nBft/a 
nB o 

rn = nB + g 

f, = sinX cos'l' - — cosX 

e e K nB r e e 

A 

and L <7 = am P litude of periodic blade loading defined in 

y Eq. (3.1) for u>. = 0 , evaluated at r = r . 

n e 

We note that the harmonic index n in the above equation 
takes on both positive l and negative integer values , and the 
absolute value of the right hand side is ‘implied. The index 
g also takes on both positive and negative integer values 


according to Eq. (2.3) describing the inflow distortion. When 
a rotor is subjected to a known inflow distortion, the re- 
sulting sound pressure can be obtained by summing over g, the 
result of Eq. (A2 . 6 ), taking into account the amplitude and 
phase angle of each gth component, before taking the rms 
value. To indicate the possibility of high levels of blade 
passing frequency noise from the small scale fan, we substi- 
tuted into Eq. (3.7) the blade loading L corresponding to 
^g }l = 0 • 01U a sinX e with U a =83.5 ft/sec , for all values of 
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the index q. and r = i-(r,+r, ). For each value of a the results 

e 2 h t a 

obtained for n = +1 and -1 from Eq. (3.7) are added to give 
the sound pressure at the blade passing frequency. The 
results of our computations for \g\ from 6 to 15 are shown 
in Fig. 3.1. For values of \g\ less than 6, the sound pres- 
sure levels computed are too low to be shown in the figure. 
Computations for |g| between 15 and 24 were not carried out, 
but an examination of Eq. (3.7) indicates that the results 
could be comparable to those presented in Fig. 3.1 when the 
correspondence of the values of m is recognized. However, as 
\g\ increases the blade loading is reduced and also the 

directivity would be affected due to the second term in f, . 

e 

It is interesting to note that even a small amplitude 
of the g’th component of inflow distortion can yield consid- 
erable radiated sound pressure, as the value of g is increased. 

I 

However, the inflow distortion over the rotor annulus has to 
be quite large to yield the amplitude of one percent of axial 
flow for the g th component used in our acoustic computations, 
as the value of g is increased. 

The special case of acoustic radiation from the rotor 
blades operating in uniform axial flow can be obtained by 


replacing in Eq. (3.7) by the steady state blade lift L q 
H ence, for the rotor in uniform axial inflow we have 


. K 

p = — 2 £_. B l *f, *J n (K r sinTMr. - r, ) 
n 4/2 ttR 0 X e nB nB e t h 


(3.8) 
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(3.9) 


where I» q - X P 0 V ^ C L c , evaluated at r = r g 

and n takes on both positive and negative values . 

( 9 ) 

Gutm , nearly four decades ago, investigated noise 
from a rotor in uniform flow and the above equation for tone 
noise is in agreement with his results. In carrying out 
acoustic computations on the small scale fan according to 
■ above Eq. (3.8), we used U =83.5 ft/sec arid C T = 0.3, in 

view of the lightly loaded condition of the rotor blades. The 
results of such computations are shown in Fig. 3.2 and these 
theoretical estimates are several orders of magnitude lower 
than test data as can be expected from multibladed rotors. I 

i 

| . 

3*2 Effect of Blade Thickness j ' > 

Thej distribution of volume sources along the chord of j 

1 I 

the jth blade to account for the blade thickness is described 
by Eq. (2.14) derived in subsection 2.3. In the presence of 
inflow distortion, the relative velocity V r met by the blade 
element, traveling along the circumference 27 rr, will be a 
function of time with its ghth component defined by Eq. (2.4). 
Consequently, we can replace Eq. (2.14) by 
db , . 

q gh,u - dT d ^ dr ’V jr ’ ex P { - i(a Vz + ignTj +1*^} ( 3 . 10 ) 

where <J>^ represents the phase angle of U k ^ and t . represents 
the time lag with respect to the "first" blade as given by 
Eq. (3.2). we note that the various functions on the right 
hand side are not identified with subscript j since all the 
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blades are considered identical, except for time delay 
accounted by x . . The location of the periodic volume source 
described above is given by Eq. (2.15) of the preceeding sec- 
tion. By using complex Fourier series representation of the 
above volume source at its coordinates, the "monopole" per 
unit volume d? dr rd9 can be obtained as 

+ im (? sinX ) /r - i(m -g) (j - 1) 2 tt/B} (3 

We note the above expression is for the contribution from 

the volume source q , „ . at £ on the chord of the jth blade, 

?«7 

due to the effect of the ghth component of the inflow dis- 
tortion. 

The far field acoustic pressure due to the elemental 
volume source, described above and located at (-? cosA,r,0) 
can be written as 

r s P 0 

KC ~ 4irS 

where the square bracket denotes that the derivative is 
evaluated at retarded time (t-S/a Q ). Substitution of Eq. 
(3.11) into the above, and application of the far field 
approximations are given in section A4 of the appendix. By 
including the volume sources distributed over the chord of 
each blade element and all the blade elements in the fan 
annulus, we obtain the far field sound pressure due to the 
ghth. component of inflow distortion as 


a 

at 


Sq 


gh , U 


(3 


/ f B C/ / 

-// 5 / 


gh y t,j 


(3.13) 


r h 0 


After integrating over 6., and summing over j, it is shown in 
section A4 of the appendix that tone noise at blade passing 
frequency exists only when = 0, i.e. for steady state 
inflow distortion. In terms of the representative radius r gJ 
the rms value of the nth harmonic of sound pressure given by 
Eq. (3.13) is derived in the appendix as 
p nBft ^ 

p _ = — -U *B*G *J (K r sin^ ) • (r . - r, ) (3 

g 4/2itR g > T mge m nB e t h 


where K nQ = nB^/a ( 

m = nB + g 
c/2 


and G 


-c/2 


• exp{im(£ sinX)/r 


+ iK nB (£ cosX ) cos'i' } *d£ ; at r 


We observe that evaluation of sound pressure according 
to the above Eq. (3.14) for the effect of blade thickness is 

complicated since the function G has to be obtained at 

mge 

each azimuth angle and harmonic index n, for a given blade 
profile. As a typical example we considered a double- 
circular arc profile with 10% thickness ratio and the evalu- 
ation of G is described in section A4 of the appendix. 
mge 

In carrying out acoustic computations on the small scale 
fan as per above Eq. (3.14), we employed r = 0.0lU a cosX e 
with U =83.5 ft/sec for all values of \g\ from 6 to 15, and 
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the resulting sound pressure levels are shown in Fig. 3.3. 
The influence of \g\ is to some extent similar to that dis- 
cussed in the preceeding subsection. We observe that the 
blade thickness effect can be comparable to that of blade 
loading, in the prediction of rotor noise due to inflow dis- 
tortion. 

For the case of the rotor in uniform axial inflow. 

Eg. (3.14) reduces to 


p nBS) 

p S = — V B* G _ *J „(K „r sint) • (r, - r. ) 

n 4/2 -itR r nBe nB nB e t h 


(3 


where G nBe = integral of Eq. (3.15) with m = nB. 

Results of computations according to above equation for the 
design operating conditions of the small scale fan are shown 
in Fig. 3.4, and we observe that these predicted values also 
are several orders of magnitude lower than test data, indi- 
cating that consideration of uniform axial inflow leads to 
insignificant noise levels. 


3.3 Quadrupole Radiation from Fluctuating Reynolds Stresses 

Our primary interest in the present investigation is on 
the noise source in a rotor that can give rise to the high 
tone levels measured at on-axis point. Consequently, we 
limited our attention only to the longitudinal axial quad- 
rupoles, by setting i = j = 1 in Eq. (2.16) for the stress 
tensor. The far field acoustic pressure from such elemental 
quadrupole can be written as 


.16) 
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• d vol 


(3.17) 


6P q = 225!ipi- (p v v ) 

4tra 2 S [_3 1 2 0 1 1 . 

o 

where the square bracket denotes that the derivative is 

evaluated at retarded time t - S/a . 

o 

Since we are interested only in the consequences of a 
^/jth component of inflow distortion, the axial component V 
in the region upstream of the rotor can be written as 


V, = U , + U (3 

1 gh } i P, i 

where U ^ i can be described in terms of its amplitude and 
phase angle in a manner similar to Eq. (2.3) or (2.4) and 
Up is described by Eq. (2.5). Considering the time deriv- 
ative in Eq. (3.17), we observe that tone noise at harmonics 


of blade passing frequency would occur only when co^ = 0 and 
we can substitute 


% V ,V 2 »o V, •V,„- exp(i( 9 e 

+ inB (0 - qt) + i(J> pi (3 

into Eq. (3.17) to obtain the nth harmonic of the blade 
passing frequency of 5p^. By including the radiation from 
all the quadrupoles present in the region upstream of the 
rotor, we obtain 

PgU = /// 6pq n rdr de dr) i (3 

r e n j 

Substitution of Eq. (3.19) into Eq. (3.17) to obtain the nth 
harmonic Sp^ and the details of integration are given in 
section A. 6 of the appendix. From Eq. (A6.2) for p q n we can 


.18) 


.19) 


. 20 ) 
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obtain its rms value as 


K ^ 

aq = , p 0, _ n _ B . cos 2 T e_ . 

9,n 4/2ttR |n| 


A 

2U 

g> i 


A 

u 

Pun 


• J (K 
m 


nB 


r 

e 


sin'}') 


X 


r e (r t 




(3.21) 


where m = nB + g 

and K _ = nBfi/a . 
nB o 

We note that the above equation is a result of considering 

only the longitudinal axial quadrupoles , and for azimuth angles 

f / 0 , one should include the effects of U and U also. 

P3., n g t 3 

Consequently any comparison of the above result with the corres- 
ponding value for dipole radiation from blade loading can be 
made only for on-axis point. 

In evaluating the sound pressure at on-axis points from 
Eqs . (3.7) and (3.21), we note that 


g = -nB 

Using the approximate expression 

|S(y) I = 


TrnBqc 


->5 


for Sears' lift response function to obtain L for substitu- 


tion into Eq. (3.7), we obtain 


{’?» < } 


= 2/2 


M. 


at ¥=0 


g 


|u u , , 

S-2-k — P . U . n .^ (sinX 

10 . V 

g» l r 


(3.22) 


The above ratio denotes the relative importance of the quad- 
rupole radiation in estimating on-axis noise from rotors in 
distorted inflow. Even though g = -nB for the ^th component 
of inflow distortion in evaluating the on-axis noise, we 
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retained both the subscripts g and n in the above equation. 

We note that all other factors except M • (U /V ) and the 

r pi,n r 

harmonic index n can be of the order one in most fan designs 
including the small scale fan. Substituting the value for 

A 

U p b n estimated on the small scale rotor in reference 3, 
we find that the ratio given by Eq . (3.22) is low. Hence, 

further consideration of quadrupole radiation from the small 
scale rotor in distorted inflow is ignored. However, for 
evaluating tone noise at higher harmonics from fans with 
blades operating at = 1.0 and M^_ = 1.0, the quadrupole 
radiation cannot be ignored. 
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4.0 


DISCRETE TONE NOISE FROM ROTOR INTERACTIONS 
WITH INFLOW TURBULENCE 

The effect of turbulence in the inflow to the rotor is to 
cause random fluctuations in the relative velocity and the 
incidence angle met by each blade element as it travels around 
the rotor circumference . Even though the events occurring on 
each blade element are random in nature, the analyses presented 
in the following two subsections show that discrete tones can 
appear in the spectrum of radiated noise from the rotor. These 
discrete tones occurring at the blade passing frequency and its 
harmonics are caused by the low frequency fluctuations of the 
turbulent inflow, which in turn are related to the longitudinal 
length scale in the frozen convected model used in our analysis. 
Because of the periodicity of the rotor-induced velocity fluctu- 
ations, quadrupole radiation from the fluctuating Reynolds 
stresses also gives rise to discrete tone noise as discussed 
in the last subsection. 

4.1 Effect of Random Blade Loading 

The fluctuating lift on a blade element is related to the 
fluctuations of u x , the turbulent velocity component normal 
to the blade chord. The latter, being a stationary random 
function of time with zero mean, can be represented by means 
of a Fourier-Stielt jes integral ^ 1 0 ^ . 

CO 

/ -iu>t 

e dZ u (to) (4.1) 

0J = -oo 
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where Z (go) is a non-dif f erentiable random function, and its 
power spectral density is given by 


dZ (go) *dZ* (go") = <£ (go) • 6 (go - go ") ‘dco dto" (4.2) 

u i u i u x 

Although u x is a random function of time and space, it has 
coherency over a certain volume in space, called the eddy 
volume. Let us assume that u : is perfectly correlated over 


a circumf erential length fc 3 . Since this velocity pattern 
repeats itself in the 9 direction with period 2 it, we can 
write for the velocity fluctuation at any 6 location due to 
the influence of one particular eddy located at 0, and radius 
r , as 


u j (r , 0 , t) = £ 

9 


A e 
i 9 


i<? (6-0 0 ) 


OQ 

/ 


-iu)t 


dZ (go) 
u x 


(4.3) 


where A = — sin(^£ 3 /2r) 


G 0=-“> 

g 9 TT a ^/2r) . 

. j j 

Since a blade element travels through angle 0 = fit 
during a time t, it experiences velocity fluctuations at cir- 
cular frequency (go - gSl) . The fluctuating lift force per unit 
span on blade element of the "first” blade located at radius 
r is then given by 


-i^9 0 r -i(Go-^fi)t 

AL 1 = 7rp o V r C 52 A g e J e £( Y')dZ u (go) (4.4) 

9 G0 = — 00 

where AL = the load due to random velocity fluctuations 

occurring in a coherent eddy, and the subscript 
1 denotes the "first" blade, 

S(y')= Sears' lift response function 

y' = reduced frequency = {w -gQ) . 
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Similar to the above Eq. (4.4), the load on the jth blade 
can be written as 


AL. 

3 


irp V c V 
o r 

9 

CO 

x / 

0J=-°° 


A *exp{-i<7(0 - Qt ) } 
g * a o j 

-i (ta-gQ) t 

e Siy'Y&Z (to) 

u. 


(4.5) 


where t. = (j - l)2iT/Bn 

3 

is the phase difference from the loading on the "first" blade. 
At any time t, the location of the point of application of 

the force L. on the jth blade element is defined by Eq. (2.11) 

3 

given in the earlier subsection 2.2. Using complex Fourier 
series representation, as employed in deriving Eq. (3.4), the 
random lift load on an elemental area rd0 dr at (-£' cosX,r,0) 
due to the loading on the jth blade, can be written as 
p V c 

6L . = 0 j — Y' A exp{im0 + i ml,' sinX/r - igO } 

m g y 


O t r 

X exp{-i (m - g) (j‘ - 1) -g-} 


00 


X 



exp{-i (tJ-g'^+mQ)t}s(Y > )dZ u (w) drd0 


0 )=-°° 

The acoustic pressure at the far field point due to the 
above impressed force can be obtained from 


6p ? = -i — i — [ -^r- { (-sinX cos'? + cosX sin0. sinf) *6L. } 

J 47rSa 0 L 9t 3 3 J 

where the square bracket indicates that the derivative is 
evaluated at retarded time t - S/a Q . Substitution of 


(4.6) 


(4.7) 
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Eq. (4.7) into the above, and far field approximation is 
given in section A3 of the appendix. After considering the 
loading on all the blade elements contained in the fan annulus, 
and the contribution from all independent coherent regions 
present in the fan annulus, the following expression is de- 
rived in the appendix for the spectrum of the far field sound 
pressure . 


4> d (v) 


-*/ 


( 7rp o V 


c) 2 -( 


V 


4iTRa 


2ur 

'T— 

3 


x £ H A a * f X * J m ^ r (u)-\S(y')\ 2 -dr (4.8) 

m g y o U J- 

where cj> (w) = u 2 (£/ttU ){1+ (JUo/U ) 2 } 1 
u^ a a 

w = v - nBf! ; m - nB + g ; n = ....-1,0,1 

S(y') = Sears' lift response function 
y' = (to - gtt) (c/2V r ) 
and f^ = sinA cos’l' - (ma Q /vr) cosA 

It is to be noted that, in the representation employed 
in Eq. (4.3) the coefficient A^ can be zero for certain values 
of £ 3 . This particular behavior of A^ is due to our assump- 
tion of a perfectly correlated eddy which of course is an 
oversimplification. Since the eddy size H is not unique, 
instead of defining A as in Eq. (4.3), it is more approp- 
riate to use the envelope function of A 2 . A reasonable approx- 
imation to this envelope function is 
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for 


A 2 

9 


sin (gSL/2r) } 2 



< 


IT 

2 


and A 2 = {(£ /2TTr ) (1 + a£ /2r)} 2 {l+(o£ /2r) 2 } 2 otherwise (4.9) 

g 3 3 3 

In carrying out acoustic computation on the small scale 
fan, the integration over r in Eq. (4.8) is approximated by 
considering the value of the integrand at an effective radius 
r g , as in the previous section 3. Assuming 3% turbulence 
intensity and parametric values of £/d e ratio, the sound spec- 
tra at an axis point 5 ft. from rotor are computed, and the 
results are given in Fig. 4.1. As the turbulence length 
scale is increased, the spectra show prominent peaks at blade 
passing frequency and its higher harmonics . We note that 
there are no spurious peaks in the spectrum obtained with the 
present approach as compared with that reported in reference 11, 
wherein perfect .correlation of loading on a finite number of 
blades was assumed. The function of 6 , the number of blades 
"chopping" an eddy, as postulated in reference 11, was respon- 

I 

sible for the series of less prominent peaks seen in the noise 
spectra. In the present approach, we considered all the rotor 
blades, but only that part of the loading caused by velocity 
fluctuations in a coherent eddy. 

Directivity patterns of radiated acoustic pressure at 
blade passing frequency and its first harmonic, calculated 
for the case of £/d e = 5.0 are presented in Fig. 4.2, and 
indicate that the noise on the side line can be 10 dB lower 
than that along the axis. 
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. . . ( 1 2 ) 

Computations carried out by Mam for acoustic power 

radiated from a rotor interacting with inflow turbulence 
indicate discrete tones at blade passing frequency, and higher 
harmonics, similar to the results presented here. However, 
the peaks in the spectra presented in fig. 4.1 are much 
sharper than the results given in reference 12, possibly due 
to the differences in the assumption of loading occurring on 

( i 3 ) 

the blades. Recently, Homicz and George presented an 

analysis of acoustic radiation from rotors in turbulent flow. 
Their results also do not show, probably for the same reasons 
as indicated above, the sharpness of the discrete tones evi- 
dent in the narrow band analysis of experimental measurements. 

4.2 Effect of Blade Thickness 

The strength of the volume sources distributed along the 

blade chord to represent the thickness effect depend upon the 

relative velocity met by the blade as discussed in the earlier 

subsection 3.2. Let us consider the consequences of the random 
i ■ , | 

nature of these volume sources resulting from the turbulent 

i , 

fluctuation in the inflow to the rotor. 1 

Similar to Eq. (4.1) we can represent u r , the fluctu- 
ating velocity component along the mean relative velocity 
direction as 
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and replace V r in Eq. (2.14) by the above expression to 
obtain the random fluctuation of the elemental volume source 

q . located at £ on the chord of the jth blade. Using the 

h» <7 

same approach as employed in the preceeding subsection, we 
obtain the randomly fluctuating volume source distribution, 
due to the velocity fluctuations in a coherent eddy, as 

6q C 7 = 2tF £ X) A *exp{im6 + im£ sinX/r - igQ } 
w m g y 

X exp{-i (m - g) (j - 1) 

OO 

X J exp{-i (oj - gQ + mft) t) d£ dr d6 -dZ u (w) (4.11) 

CJ=-00 r 

Acoustic radiation from above random volume source 
distributuion can be calculated from 


6p s . 

P £J 


4irS 


3T ^ 


(4.12) 


where the square bracket is used in the same sense as in 
Eq. (4.7). The substitution of Eq. (4.11) into Eq. (4.12) is 
shown in section A5 of the appendix and the acoustic spectrum 
due to radiation from all the blades in the fan annulus is given 

b Y r 

, r if v p 1 2 

14# -V(2*r/£ 3 ).B 2 

r. 


X E E ^-|G 

mg y 


mge^^u (U),J ^ ( ^ r sin,,)dr 


(4.13) 
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where <$ (to) = spectrum of u as defined by Eq. (2.8) 
u r r 

to = v - nBft 

m = nB + g 


and 


^ - 
9 


= function defined in Eq. (4.9) 


mge 


c/2 

7 # 


-c/2 


)exp{im£ sinX/r + i — £ cosX cos^l-d? 

cl 

0 


evaluated at r = r 


(4.14) 


In carrying out acoustic computations on the small scale 


fan as per above Eq. (4.13), we employed the same turbulence 
parameters as used in subsection 4.1. The integral in 
Eq. (4.14) is evaluated for the same blade profile used in 
the computations presented in subsection 3.2. The computed 
spectra presented in Fig. 4.3 also show prominent peaks at 
blade passing frequency and its higher harmonics but the 
levels are much lower than those given in Fig. 4.1, for the 
blade loading effects. The directivity patterns of tone 
noise at the blade passing frequency and its first harmonic, 
as computed from Eq. (4.13) are shown in Fig. 4.4. 

Computational results presented in this subsection 
indicate that blade thickness effect can be important at 
higher harmonics of the blade passing frequency. 


4.3 Radiation from Random Quadrupole Sources 

In the region upstream of the rotor leading edge there 
exist periodic velocity fluctuations, as discussed in sub- 
section 2.1.2, in addition to the random turbulent velocity 
fluctuations of the inflow described in subsection 2.1.3. 
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Since our primary interest is to evaluate such sources that 
would yield acoustic radiation to on-axis points, we shall 
consider only the longitudinal axial quadrupoles as in sub- 
section 3.3. For random quadrupole calculations, instead of 
using Eq. (3.17), it is more appropriate to obtain the spec- 
trum of the far field acoustic pressure by integrating the 
space-time correlation of the fluctuating Reynolds stresses 
over the flow region, as indicated in reference 14. Thus, 
the far field -sound spectrum due to the longitudinal axial 
quadrupoles can be written as 

<$^(R,v) = p2 o cos4y y 1 * f d 3 n /d 3 c 

(4Tra 2 R) 2 2 tt ■' •' 

o OO oo 


/' 


R(n,5,T) • COS (T. - T*)dT 


(4.15) 


where R(n,5,x) = V 2 (n , t) V 2 (rj + § , t + x) (4.16) 

Vj = axial component of V described in Eq. (2.1) 

r 

and t* - £*R/a Q R is the retarded time difference. 

The position vector n , the separation vector £, and the region 
of integration are schematically 1 shown in Fig. 4.5. ’In con- 
sidering the tone noise of the quadrupole radiation due to 
inflow turbulence and rotor induced potential velocity fluc- 
tuations, we need to retain only the product term 2U u of V 2 

pi i i 

in Eq. (4.16). The derivation of the space-time correlation 
function and the limits of integration are discussed in 
detail in section A7 of the appendix. After certain amount of 
algebraic simplification, the spectrum of far field acoustic 

pressure due to radiation from quadrupole sources can be 
written as 


36 


. P v‘ 

_,v) ° 


pn 


|47ra 2 R. 


0 0 

where 

A = fan annulus area 


12 A(<p - $ ) £ 4 • (u 2 ) • (u ) 2 
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■nl I (U k ) 2 + (w - V) 2 (U k ) 2 + (w + v) 2 
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k = n2iT/d , on = nBQ. 

P P 


For the same turbulence parameters, employed in calcu- 
lating the results given in Fig. 4.1, we evaluated the quad- 
rupole radiation to the on-axis point from the above equation, 
using the fundamental of the rotor-related potential velocity 
fluctuations. The resulting spectra presented in Fig. 4.6 
peak at blade passing frequency , and once again the sound 
pressure levels are far below the corresponding values shown 

I 

m Fig. 4.1 i for radiation from blade loading. The spectra 
shown in Fig. 4.6 peak only at the blade passing frequency, 

t 

. • j 

since the fundamental component of U was considered. Based 

Pi 

on his computations for tone power estimates at blade passing 
frequency, Chandrashekhara ^ 1 5 ^ also reached similar conclu- 
sions. 

We observe that on the right hand side of Eq. (4.17), 
the dominant influence is from the first term of , which 
c an be recognized as the inflow turbulence spectrum shifted to 
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the frequency a> . This feature of quadrupole radiation due to 

( 16 ) 

inflow turbulence was also noted by Williams and Hawkings 
The various other terms containing k^ , not discussed in ref- 
erence 16 , appear in the above equation due to the exponential 

/s 

decay of U with distance upstream of rotor and integration 
I? 1 

over the whole flow region. A comparison of the Eqs . (4.17) 

and (4.8) evaluated at on-axis points can be reduced to 


4M 2 sinX (u 


at Y = 0 


Since U can never be larger than V , the on-axis tone 

pi,n 3 r 

noise from subsonic rotors in turbulent inflow is dominated 


by the dipole radiation at all harmonics. 
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5.0 EXPERIMENTAL MEASUREMENTS AND COMPARISON 

WITH THEORETICAL PREDICTIONS 

In the experiments conducted by NASA, Ames Research 
Center, the small scale rotor, whose design parameters were 
employed in the acoustic computations presented in the pre- 
ceeding sections, is located at the center of the anechoic 
chamber and the^exit flow from, the fan is vented out of the 
chamber to avoid recirculation into the rotor inlet. Flow 
measurements , taken with the hot-wire located in the fan 
entrance bell-mouth at the pitch radius, and 6 in. upstream 
of the rotor, indicate a mean axial flow of 83.5 ft/sec and 
turbulence intensity of 0.75%. The spectrum of the axial 
component of turbulence, obtained with a 50 Hz band-width, 
is shown in Fig. 5.1. From the autocorrelation function 
presented in Fig. 5.2, a longitudinal velocity correlation 
length scale of 20 inches is deduced. 

In carrying out acoustic computations consistent with 
the test conditions we assumed "f rozen-convected" turbulence 
hypothesis, with 83.5 ft/sec as the convection velocity to 
obtain the power spectrum of u from Eq. (2.7). Based on 
the measurements reported by Chandrashekhara ^ 1 5 ^ on a small 
scale rotor of comparable size, we find that the length scale 
H in the circumferential direction can be much smaller than 

3 

20 inches. Consequently, we employed parametric values of 
£ 3 = 1, 2 and 3 inches along with the above mentioned spec- 
trum for u : to compute the noise spectra, according to 
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Eq. (4.8), at 7 ft. distance along the axis. The results thus 
obtained for dipole radiation from blade loading, are shown in 
Fig. 5.3. The peak levels at blade passing frequency and its 
harmonic did not vary appreciably with Z 3 in the range examined. 
Similar results obtained from Eq. (4.13) for monopole radiation 
from blade thickness effect are shown in Fig. 5.4, indicating 
the minor role of blade thickness in rotor noise prediction at 
on-axis points. Based on the conclusions reached in preceeding 
section regarding the relative unimportance of noise from the 
quadrupole sources, computations using Eq. (4.17) were not 
carried out for the test configuration. 

The noise measured at 7 ft. distance from the rotor, along 
its axis, is analyzed using a 50 Hz bandwidth. The resulting 
spectrum, shown in Fig. 5.5, contains sharp peaks at the blade 
passing frequency and its first harmonic. For correlation with 
theoretical estimates , we note that the latter are computed 
from the sound pressure spectrum given by Eq. (4.8), which is 
derived for both positive and negative values of radian fre- 
quency. Hence, the acoustic pressure spectral density per Hertz 
for positive frequencies is obtained by using a multiplying 
factor of 4 tt to the results of Eq. (4.8). The spectral density 
thus obtained is integrated over 50 Hz bandwidth before finding 
the dB level for comparison with the measured values. Using 
length scale Z 3 = 3 inches, we carried out such calculations in 
the neighborhood of the blade passing frequency and its first 
harmonic and the theoretical estimates are shown by dotted lines 
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in Fig. 5.5. The change in the scale of the abscissa and inte- 
gration over the 50 Hz bandwidth are responsible for the dif- 
ference in estimated sound pressure levels presented in Fig. 5.5 
and those shown in Fig.. 5.3. It appears that the sharpness of 
the discrete tones observed in the test data is predictable 
from the theory presented here. Variability of the measured 
tone levels obtained with a 50 Hz bandwidth filter centered at 
—■ -- -the blade passing frequency and its -first harmonic are shown 

in Fig. 5.6. Our theoretical predictions of the tone levels 
indicated by the dotted lines in Fig. 5.6 appear to be in 
reasonable agreement with the measured values. 

The variability of the tone level, with time as seen in 
Fig. 5.6, can have its origin in the very nature of the noise 
sources. An examination of Eq. (4.8) with ¥ = 0 , indicates 
that the tone noise level is proportional to the mean square 
value of turbulent velocity fluctuations at the low end of the 
spectrum. On the other hand, the data presented in Fig. 5.6 is 
obtained using a short averaging time-constant in view of the 
high frequencies of the tones. Consequently, the low frequency 
nature of the generation mechanism appears in the output as 
variability of the tone level. 

The inlet to the rotor is not yet surveyed to assess the 
steady state distortion of the inflow. Hence, no estimates 
can be made at the present time for the on-axis noise from 
the ordered sources discussed in section 3. 
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6.0 


CONCLUSIONS AND RECOMMENDATIONS 


The theoretical investigations of noise sources in a 
rotor and typical computational results presented in this 
report lead us to the following observations: 

(a) Steady state distortion of inflow, even of small mag- 
nitude, can generate high levels of noise at blade passing 
frequency harmonics. 

(b) Periodic inflow distortions lead to tones at frequen- 
cies other than the blade passing frequency harmonics. 

(c) The interaction of the rotor blades with inflow tur- 
bulence gives rise to noise whose spectrum contains prom- 
inent peaks at blade passing frequency and its higher 
harmonics. The sharpness and level of these discrete 
tones depend upon the intensity of turbulence and its 
velocity correlation length scale. 

(d) Computations presented in this report indicate that 
the discrete tone level increases as the ratio i/d is 
increased. Since blade speed and the axial velocity 

U through the rotor are held constant in our computations, 

a 

increased value of i/d can be interpreted as more blades 
"chopping" a coherent eddy as it is convected across the 
rotor. 

(e) Theoretical investigations show that the origin of 
the discrete tone noise lies in the turbulent velocity 
fluctuations at the low-end of frequency spectrum. 
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(f) Limited computations for the effect of inflow tur- 
bulence on the small scale rotor indicate that blade 
loading fluctuations are the dominant sources of noise in 
isolated rotors. 

It is necessary to include the acoustic radiation from 
all the quadrupole sources and also the influence of the duct to 
make a rational comparison of the various noise sources in 
a rotor. 

The analyses of noise sources and radiation therefrom, 
presented in this report, are applicable to rotor blades 
operating at relative Mach numbers in the low subsonic range. 

To evaluate noise from lift fans of present day configurations, 
the analyses must be extended to the consideration of high 
relative flow Mach numbers. 

Space-time correlation measurements of inflow turbulence 
for various inlet configurations are recommended to determine 
the characteristics of turbulent velocity fluctuations. Such 
flow measurements along with corresponding acoustic data 
would confirm the assumptions employed in the present analysis 
and provide a basis for further improvements in the theoretical 
prediction of rotor noise. 

From the preliminary investigation of noise sources 
presented here it appears that reduction of discrete tone noise 
from rotors can be achieved by eliminating inflow distortion 
and decreasing the length scale of the turbulent velocity 
fluctuations . 
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APPENDIX 


The steps involved in deriving the expressions for the 
far field acoustic radiation from a rotor subjected to vari- 
ous inflow conditions are given in detail in the following 
sections of this appendix. 

Al Coordinate System and Distance to Far Field 

Let R and S denote the vector distance of the far field 
point from the origin and source respectively as shown in 
Fig. 2.1. For the sake of convenience and considering the cyl- 
indrical symmetry of the radiated field, we have chosen the 
meridional plane 0 = 0 to pass through the field point. By 
appropriate translation and rotation of the coordinate system, 
we can obtain the components of S in the axial, radial and cir- 
cumferential directions as 

S = R cosf -x 

1 

S = R sin¥ cos0 - r 

2 

S = -R sin¥ sin0 (Al.l) 

3 

The distance of the field point from the source location can 
be approximated as 

S = R - r siny cos0 - x cos'F (Al. 2) 

In the far field computations of fan noise the second order 
terms r/R and x/R can be ignored in the amplitude decay term 
but not in the retarded time function. 
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A2 


Dipole Radiation from Periodic Blade Loading 


Substituting Eq. (3.4) into Eq. (3.5) and using Eq. 

(A1.2) for the far field approximation with the corres- 
ponding values of x. and 0 . from Eq. (3.3), we obtain the 

0 0 

following equation for acoustic radiation from an elemental 
area rd0 dr of the rotor disc 

, L , *drd0 iK , 

8p , . = -& - = ■ (sinl cos'F - cosX sine sinV) 

2ir ^ 4ttR 

X exp{-iK ,a (t - R/a ) + i<|> ,} 
c mgh o o gh 

X exp{im0 - iK ,r sin^ cos0} 

mgh 

X exp{im U' h sinX)/r + iK^ cos X) cosy} 

X exp{-i (m - g) (j - 1) 2 tt/B} (A2.1) 


where K . = (w, - + mfl) /a 

and X = blade stagger, employed in obtaining axial and 

circumferential components of force L^. 

On the left hand side of the above equation the subscript gh,j 

indicates that it is the contribution from the loading on 

the jth blade element due to its interaction with the ghth 

component of the inflow distortion. Considering the loading 

on all the blade elements contained in the fan annulus, we 

obtain the far field acoustic pressure as 


2t r 


d 

p gh 


■f fk 


5p 




(A2.2) 


'h 
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The superscript d and the subscript gh on the left hand side 
denotes that the expression is for the contribution from dipole 
radiation from the effect of only the ghth component of the in- 
flow distortion. Substituting Eq.(A2.1) into Eq.(A2.2) and 
using the following identities 
2tt 


/ 


2tt 

J C sin8 *exp{i (m8 - zcos8 )}d0 = {—■) 2tt (-i) m J m (z) and (A 2. 4) 
o 

B 

£ exp{-i {m g) {j - 1) 27r/B} = B j if = an integer 

J=1 

= 0 otherwise (A2.5) 

we obtain the nth component of far field sound pressure as 
r , 


exp{ i (m0 - zcos0 )}d8 = 2t\ {-\) m J (z) 

171 


(A2.3) 


^ <- i >" 


iK ,a (t - R/a ) + i<j> ,}*f. 

o o gn X 


X exp{im(Zg h sinX)/r + iK cosX)cosf} 

x I V J „ ( V- r sin,,, - ar 

“ here K m s7 f. - 


f . . = sinX cos’i' - ^ 

A i\ 


77 ? 


‘mgh' r 


co sX 


and 


m = nB + g 


(A2 . 6 ) 


We note that the frequency of the nth component of acoustic 

pressure defined above is 

. (.m - g)tt 

— — = — + nBN , 

2 TT 2tt 


(A2.7) 


where N = rotor rps . 
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Consequently, sound radiation would occur at blade passing 
frequency only when = 0, i.e. for steady state inflow dis- 
tortion. By using a representative radius r g and replacing 
the integrand of the above Eq. (A2.6) by its value at r g , 
we can then write the rms value of the nth harmonic as 


K 


nB 


g n 


where K 


nB 


-4/2ttR 

nB£2/a 


BL • f. • J (K _ r sinf) *(r. -r, ) 
g Xe m nB e t h 


(A2 . 8 ) 


and L 


amplitude of periodic blade loading defined 
by Eq. (3.1) for = 0 . 


A3 


Dipole Radiation from Random Blade Loading 


Substituting Eq. (4.6) into Eq. (4.7) and using 
Eq. (A1.2) for the far field approximation we obtain the 
following expression 

, p V c r iK 

6p . = - ° - y — 5Z H ^ s; '- n ' >k cos'? - cosX sin0 sin’!') 

^ m q & J 77 

» u =-oo 

X exp(-iK a (t - R/a ) - i gQ } 
c mg o o v o 

X exp{im0 - iK r sin¥ cos0} 
mg 

X exp{i m£' sinX/r + iK t,' cosX cos^} 

mg 

X exp{-i (j - 1) (m - g) 2ir/B} S(Y")dZ (w) dr d0 (A3.1) 

u i 

where K = (w-gSi+mtt )/ a 
mg o 

S(y') = Sears' lift response function at reduced fre- 
quency y , 

and y' = (w - gQ) ■£- 
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Now carrying out summation over j from 1 to B to account for 
the effect of all blades and integrating over 0 <_ 6 <_ 2tt 
gives the far field acoustic pressure from blade elements 
within an annulus 2irrdr due to interactions with a coherent 
eddy as 


p V cB 
, d r 

Ap = 


oo 

£ £ / 

mg * J 


iK 

""4rrR"" 2Tr (~i) m 


0)=-a> 


X exp{-iK a (t - R/a ) - i gQ } 
mg o o a o 

X exp{im?' + iK r' cosX cosy} 

r mg 


X J„,(K_ r siny) -f x *5(Y^)dZ u ^ (co)dr 


m mg 
where m - nB + g , 


and f x = sinX cosy - 


cosX 


(A3. 2) 


mg 

Noting that the dimension of the coherent eddy in the radial 
direction is l 2 we can approximate the integration with 
respect to r by some mean value of Ap^ multiplied by Z 2 to 
obtain the far field acoustic pressure from blade elements 
within an annulus 27rr£ 2 . It is to be noted that the value £ 2 
in the above equation must be limited to (r fc - r^) , since the 
loading on the blade elements is defined only within the fan 
annulus. The auto-correlation function of the acoustic 
pressure, thus obtained for interactions with a coherent 
eddy, leads to the spectral density by taking its Fourier 
transform. However, it is only the contribution from the 
interaction of the rotor blades with a coherent eddy. To 
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obtain the spectrum due to all coherent regions present in 
an annulus 2irrdr, we can multiply the spectrum obtained as 
above by 2irr/Z 3 and dr/Z 2 . The spectral density at circular 
frequency v , of the radiated acoustic pressure from all the 
blade elements within the rotor annulus is then given by 


$ d (v) 


= §/<" 


o V c) 2 • ( - 

p o r c; '4irRa 


2nr 
2 i 


*^r— *B 2 


X Y* A 2 *f 2 • J 2 (— r sin! 
A g \ m a 

m g a o 


X $ u (o>) • |5(y") | 2 *dr 


(A3. 3) 


where u = v - + gQ. ; 

and y' 


m 


nB + g ; n = . . .-1,0 ,1. 


(tn-gSl) and other parameters are as defined 

r in Eqs . (A3.1) and (A3. 2). 


A4 Monopole Radiation Due to Inflow Distortion 

The effect of the g-Tzth component of the inflow on the 
thickness of the jth blade rotating with a circular frequency 
is represented by fluctuating volume sources 6q^ 
as given in Eq. (3.11). We note that the elemental volume 
rd9 dr d? referred to in the equation is located, in the annulus 
2irrdr,at coordinates given by Eq. (2.15). Substituting 
Eq. (3.11) into Eq. (3.12) and using the far field approxima- 
tions as before, we obtain the acoustic pressure due to the 
elemental monopole considered as 
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^gKU = ^ dr d9 


£ 

m 


, - lp o a A)> , 

v 4 ttR ' 


X exp{-iK^ a o (t-R/a o ) + i^} 

X exp{im0 - iK ^ r sin'i' cos8} 

X exp{im(£ sinA)/r + iK^^ (£ cosA)cos’l'} 

X exp{-i (m - g) (j - 1) 2 tt/b} 

where K . = (to, - gQ + mtt) /a„ 
mgn n o 


(A4.1) 


Including the source-sink distribution over the chord of 
each blade element and all the blade elements in the fan 
annulus, the far field sound pressure due to the ghth com- 
ponent of inflow distortion can be obtained from the follow- 
ing equation. 



6p 


gh, zj 


(A4.2) 


Substituting Eq. (A4.1) into Eq . (A4 . 2 ) and using identities 

(A2.3) and (A2.5), we obtain the nth component of sound > 


pressure as 
r t 


p S , = ( U , _ — kv mqh (~i ) m exp {-iK a (t - R/a ) + i <l> ,> 

^ gh,n J gh, r 4 ttR r mgh o o gh 


X B *G *J (K , r sin S') *dr ^ (A4.3) 

mgh m mgh 

where K „ = (co - #0+ mft)/a 
mgh h w o 

m = nB + g. 
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and 


c/2 


J mgh 


j' g^-exp{i m (c sinX)/r + iK^ cosX) cosy} -d? (A 4 . 4 ) 


-c/2 


We note that tone noise at blade passing frequency exists 
only when oj^ = 0, i.e. for steady state inflow distortion. 

Using the value at representative radius r g , for the integrand 
in Eq. (A4.3), we can write the following expression for the 
ms value of the nth harmonic of acoustic pressure due to <yth 
component of steady state inflow distortion. 


p nBft „ 

5 s = *U • B - G *J (K r sin^) • (r . - r. ) 

g n 4/2ttr 9 r mg e m nB e t h 


(A4 . 5 ) 


where K nfi = nBf2/a Q 


and 

G mg e = f gjrexpfimCs sinX)/r +iK nB (£ cosX ) cosy} -d£ (A4 . 6 ) 

-c/2 evaluated at r = r 


m = nB + g } 
c/2 


We note that the far field sound pressure given by 

Eq. (A4.5) depends upon evaluating the integral of Eq. (A4.6) 

for the specific blade section profile, at each value of 

azimuth angle and the particular values of the indicies n and 

q . To illustrate the evaluation of the function G„, , let us 

a mge 

consider the double circular arc profile described by 
db 


at “ - 4 <c/= e ) ( b /=)„, ax ,e 


(A4.7) 


where c is the blade chord and (b/c) is the thickness 

e max, e 
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ratio of the blade section, both considered at the repre- 
sentative radius r g . Using Eq. (A4.7) and the, following 
notations 


ch ' "■ 

and y = m ( j sinX) /r e + K nB (j cosX)cos¥ 

we can rewrite Eq. (A4.6) as 

1 


(A4.8) 


G = -b 
mge max 


n-ei 


yn 


dn 


(A4 . 9 ) 


From the value of the above definite integral we can show 
that 


= 2b 


l cosy siny 


1 mge 1 max ,e 1 y y 2 

where y is defined by Eq. (A4.8). 

In evaluating sound pressure at blade passing frequency at 

r 

on-axis points on the small scale fan, we note that 


(A4.10) 


leading to 


y = .28 


G = .1656 b 
1 mge 1 , max,e 


(A4.ll) 


A5 


Monopole Radiation from Random Sources 


Substituting Eq. (4.11) into Eq. (4.12) and using 
Eq. (A1.2) for the far field approximations, we obtain 
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iK a 

mg o 

47rR 


X exp{-iK m? a o (t-R/a o ) -igB Q } 

X exp{imQ - iK r sin'F cos6} 

mg 

X exptimC sinX/r + iK x, cosX cos'i'} 

mg 

X exp{-i (m - g) (j - 1) 2 tt/b} • dr d6-dz ( w ) (A5.1) 

^ u r 

where K = (ai - gfl + m£2) /a . 

■ ■ mg • " - 0 - — ~ — — -- 

Following the procedure employed in the previous subsection 
A3 step by step, we can easily derive the following expres- 


sion for spectral density of far field acoustic pressure due 
to blade thickness effect 


$ S (v) 


V 


{vp / 4 ttR) 2 SL (2nr/Z ) 

0 2 3 


B 2 


x E E Binfldr (AS. 2 ) 

mg ° r o 

where m = nB + g 

0) = V - mQ, + gtt 


and 

G 

mge 


c/2 

f if exp{im (? 
-c/2 


sinX)/r + i i \ cosX ) cost} *d£ 

a 

0 

evaluated at r = r 


e 

We note that Eq. (A5.3) is similar to Eq. (A4.6). 

Thus the amplitude of G in Eq. (A5.2) can be evaluated 

mge 

from Eq. (A4.10) and Eq. (A4.8) with replaced by v/a Q . 


(A5.3) 
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A6 RADIATION FROM ORDERED QUADRUPOLES 

Confining our discussion to longitudinal axial quad- 
rupoles only, the "fluctuating Reynolds stresses" due to 
the interaction of steady state inflow distortion and the 
potential velocity perturbation is given by Eq. (3.19). 

The far field acoustic pressure radiated by these "fluctu- 
ating Reynolds stresses" can be obtained by substituting 
Eq. (3.19) into Eq. (3.17) to give 

2 p 

6p^ = — -. Q -- nB cos 2 '}' U U exp{-2irn, |n|/d} 

^ri 4 ttR g tl p i , n 11 1 

X exp{-iK nB a Q (t - R/a Q ) } 

X exp{im0 - iK nB r sin? cose} 

X exp{ i<b „ - iK p cosTjrdr d0 dp 
Pi, n na i x 

where m = nB + g 

and K _ = nBft/a . 
nB o 

To include the radiation from all the quadrupoles 
present in the region upstream of the rotor, we can integ- 
rate Eq. (A6.1) over the region 

o < n < 00 

i 

o £ e <_ 2 7T 

r h ^ r i r t 

to arrive at the following expression for the nth harmonic 
of the far field sound pressure due to the g th component of 
the inflow distortion. 
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(A6.1) 


(d/n) 



2p K n 
o nB 

4nR 


cos 2 ’}' 



(-i) 


m 


U 

g> i 


u 

Pi,n 


X exp{-iK nB a o (t-R/a o ) + i* pi3n } 

X J (K _ r sin¥)r-dr 
m nB 


(A6 . 2 ) 


We have ignored the retarded time in the axial direction 
when performing the integration with respect to n because 
the potential velocity perturbation decay very fast in the 
axial direction. 

Using a representative radius r g and replacing the inte- 
grand of the above Eq. (A6.2) by its value at r g , we can 
write the rms value of the nth harmonic as 


p^ 

g,n 


p K „ 
o nB 

4/2?rR 


cos 


In I 




2U 


U 

g> i Pi 


•J (K .rsinflT (r. - r, ) 
, n m nB e t h 


(A6.3) 


where m = nB + g 

and K „ = nBf2/a • 
nB o 


A7 RADIATION FROM RANDOM QUADRUPOLE SOURCES 

Retaining only the product term 2U^ u in v| , we note 
that Eq. (4.16) reduces to 

^(ruS/T) = 4R U u (n,C,T) (A7.1) 

i pi i 
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Since the periodic fluctuations are not related to the 

inflow turbulence fluctuations u , the above space-time 
correlation function can be written as 


^ u = R., • R 

Pi i X. u . 


(A7 . 2 ) 


Pi i 

From the description of the nth component of U given in 
Eq. (2.5) we obtain 

It = y{U } 2 cos (k E - a) t) *exp{-2k n - k (A7.3) 

pi^n 2 P^ n P 3 P P i P 1 

where k = 2-irn/d 
P 

and a) = nBfJ 
P 


Substituting the above Eq. (A7.3), Eq. (2.9) and Eq. (A7.2) 
into Eq. (A7.1), the required space-time correlation function 
can be written as 

R v ?(n,? f T) = 2{ u pi n } 2 *u 2 «cos (k p £ - w p i) 

exp{-2k p n -k p C - |5 

U-.l/t,} »7-4) 


The above equation is then substituted into Eq. (4.15) and 
integration over the flow region is performed to obtain 
Eq. (4.17). In carrying out the integration the retarded 
time difference t* is ignored since the fluctuations U , „ 
decay rapidly with n ^ . The integration being only over 
the upstream region of the rotor, the ranges of n and £ are 
restricted by the following limits 

-n < E < 03 and 0 < n < 00 

i — x — i - 


56 


REFERENCES 


1. Rao, G.V.R. and Associates, "Study of Non-radial Stators 
for Noise Reduction", NASA CR-1882, September, 1971. 

2. Kazin, S.B. and Volk, L. J. , "LF336 Lift Fan Modification 
and Acoustic Test Program", NASA CR-1934, December, 1971. 

3. Rao , G . V . R . , et al, "Analytical Lift Fan Noise-Study", 
NASA CR-114576 , July 13, 1972. 

4. Lighthill, M.J., "Sound Generated Aerodynamically " 

The Baker ian Lecture, Proc. Roy. Soc. (London), Vol. 276, 
Series A, 1962, pp 147-186. 

5. Hanson, D.B. , "A Unified Analysis of Fan Stator Noise", 
presented at 84th meeting of Acous. Soc. Am., Miami 
Beach, Fla., Nov. 28 -Dec. 1, 19 72 

6. Deming, A.F., "Noise from Propellers with Symmetrical 
Sections at Zero Blade Angle, II", NACA TN 679, 1938. 

7. Lighthill, M.J., "On Sound Generated Aerodynamically " , 
Proc. Roy. Soc. (London), A211, 1952, pp 564-687. 

8. Kemp, N.H. and Sears, W.R., "The Unsteady Forces Due 
to Viscous Wakes in Turbomachines", J. Aerospace Sci., 
Vol. 19, 1952, pp 713. 


57 



9. Gutin, L., "On the Sound Field of a Rotating Propeller", 
NACA TM 1195, 1948 (Physikalische Zeitschrift der Sow- 
jetunion, Band 9, Heft 1, 1936, pp 57-71). 

10. Sveshnikov, A. A. , Applied Methods of the Theory of Random 
Functions , Pergamon Press, New York, 1966. 

11. Rao, G.V.R. et al, "Rotor Noise Due to Inflow Turbulence", 
A. I. A. A. paper No. 73-632, 6th Fluid and Plasma Dynamics 
Conference, Palm Springs, Calif. July 1973. 

12. Mani, R. , "Noise Due to Interaction of Inlet Turbulence 
with Isolated Stators and Rotors", J. Sound Vib. , Vol. 17, 
No. 2, 1971, pp 123-132. 

13. Homicz, G.F. and George, A.R., "Broadband Radiation from 
Rotors in Turbulent Flow", Proc. Interagency Symposium 
on Univ. Research in Transportation Noise, Stanford 
University, March 1973. 

14. Chu, W. T. , "Turbulence Measurements Relevant to Jet 
Noise", UTIA Rpt. 119, Nov. 1966. 

15. Chandrashekhara , N. , "Sound Radiation from Inflow Tur- 
bulence in Axial Flow Fans", J. Sound Vib., Vol. 19, 

No. 2, 1971. 

16. Ffowcs Williams, J.E. and Hawkings, D.L., "Theory Relating 
to the Noise of Rotating Machinery", J. Sound Vib., Vol. 
10, No. 10, 1969. 

58 





O test data 



Fig. 3.1 Effect of gth Component of Inflow Distortion on Blade Passing Frequency- 

Noise at 5 ft. from the Small Scale Rotor 
(computed from periodic blade loading) 
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Fig. 3.2 Blade Passing Frequency Noise at 5 ft. from the Small 

Scale Rotor in Uniform Flow 
(computed from steady state blade loading) 
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Fig. 4.2 Directivity of Discrete Tone Noise at 5 ft. from Small Scale Rotor 
(computed from random blade loading; 3% turbulence intensity) 
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Fig. 4.4 Directivity of Discrete Tone Noise at 5 ft. from Small Scale Rotor 
(computed from blade thickness effect; 3% turbulence intensity) 
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Fig. 5.1 Spectrum of the Axial Component of 

(band width 50 Hz) 
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Fig. 5.4 Noise Spectrum at On-axis Point 7 ft. 

from Small Scale Fan Computed from 
Blade Thickness Effect Corresponding 
to Measured Turbulence Parameters 
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Fig. 5.5 Narrow Band Analysis of Noise Measured at the On-Axis Point 

(band width 50 Hz) 
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